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Abstract 

An ideal conductor electrode in contact with a semi-infinite two-dimensional two- 
component plasma in an external potential is considered. The model is mapped 
onto an integrable sine-Gordon theory with Dirichlet boundary conditions. The 
information gained from the mapping provides an explicit form of the surface tension 
in the plasma-stability regime. 
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1 Introduction 

In a previous paper hereafter referred to as I, the bulk thermodynamic properties 
(free energy, specific heat, etc..) of a model, the two-dimensional (2D) two-component 
plasma (TCP), or Coulomb gas, have been obtained, exactly. In the present paper, a 
surface property of the same model is considered: the surface tension, at a rectilinear 
interface between an ideal conductor and the Coulomb gas, is obtained, exactly, as a 
function of the bulk density, the applied electric potential, and the temperature, in the 
whole temperature range for which the point-particle model is stable. Like in I, a mapping 
onto a sine-Gordon field theory, now with a Dirichlet boundary condition, is made, and 
known results about that field theory are used. The resulting surface tension is checked 
on its high-temperature expansion derived from a renormalized Mayer expansion and on 
its singular behaviour close to the collapse point. 

The model under consideration mimics the interface between an electrolyte (the two- 
component plasma, made of two species of point-particles, of opposite charges ±1) and 
an electrode (the ideal conductor). Classical equilibrium statistical mechanics is used. In 
the grand-canonical formalism, the control parameters are the inverse temperature /3 and 
the two fugacities z + and z_ of the positive and negative particles, respectively. Instead 
of z + and Z-, it is convenient to use z and <p defined by z± = ze ±t3ip . Alternatively, 
chemical potentials /i + and /x_ can be defined by z± = exp(/3/i ± )/A 2 , where A is the 
de Broglie thermal wavelength. The bulk properties depend only [0 on the chemical 
potential combination /i = (/i + + /i_)/2, i.e. on z, while /i + — /i_ (or ip) is relevant only 
for the surface properties ||, 0, ||. The parameter <p has a physical meaning: it is the 
electric-potential difference between the bulk and the electrode. Indeed, if the potential 
of the electrode is taken as the zero and <p is the potential in the bulk, each chemical 
potential, i.e the reversible work for adding a positive or negative particle into the bulk, 
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has an electric part ip or —ip, respectively. 

The point-particle model is stable against collapse of positive- negative pairs for (3 < 2. 
This is also the stability range in presence of a rectilinear ideal conductor wall. Indeed, 
a particle at a distance x from the wall interacts with its own image through a potential 
(1/2) ln(2x) and the corresponding Boltzmann factor (2x)~P/ 2 is integrable at small x if 
and only if (3 < 2 (at large x, the interaction is screened). 

The paper is organized as follows. In Section 2, the model is precisely defined and its 
mapping onto a sine-Gordon field theory is made. This field theory is described in Section 
3. The desired surface tension is derived in Section 4. Its high temperature expansion is 
checked in Section 5. Its singular behaviour close to the collapse point (3 = 2 is checked 
in Section 6. 

2 Mapping 

We consider an infinite 2D space of points r E R 2 defined by Cartesian coordinates (x, y). 
The model electrode-electrolyte interface is localized along the y axis, namely at {r = 
(0, y)}. The half-space x < is assumed to be occupied by an ideal conductor of dielectric 
constant e — > oo, impenetrable to particles. The electrolyte in the complementary half- 
space x > is modeled by the classical 2D TCP of point particles {j} of charge {qj = ±1}, 
immersed in a homogeneous medium of dielectric constant = 1. The interface is kept at 
zero potential while one assumes a given potential ip in the bulk. Equivalently, there is a 
splitting of the fugacities of the ± particles: 

z + = zexp((3<p), = zexp(— flip) (1) 

The system is translationally invariant in the y direction, so the position-dependent par- 
ticle densities n±(r) depend only on x. Let us denote their asymptotical x — > oo values 
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by n + = n_ = n/2 where n is the total particle number density. In the case ip = 0, 
n+(x) = n,-(x) everywhere. 

In infinite space, the Coulomb potential v at spatial position r, induced by a unit 
charge at the origin, is given by the 2D Poisson equation 

Av(r) = -2n5(r) (2) 

The solution of (Q) reads 

v(r) = -\n(\r\/r ) (3) 

where the length constant ro, which fixes the zero point of energy, is set for simplicity to 
unity. Here, the interaction energy E of particles {qj, Vj = (xj > 0, yj)} consists of two 
parts (see, e.g., [|J): 

(i) direct particle-particle interactions, 

^W(ln-^-l) (4a) 

i<j 

(ii) interactions of particles with the images of other particles and with their self-images 
due to the presence of the conducting wall, 

-^T,m<\n-r*\) (4b) 



where r* = (—x,y). Introducing the microscopic charge density p(r) = qjS(r — r 



3)i 



the energy contributions ( 4aj) and ( ^q ) are expressible as follows 

E = ^Jd 2 r J dVp(r) [v{r, r') - v{r*, r')] p{r) - ^Nv{0) (5) 
where v(0) is the self-energy. Introducing the microscopic charge + image charge density 

p( r ) = X! ij l 6 ( x - x i) - s ( x + x j)] 6 (y ~ Vi) ( 6 ) 
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the energy (|5]) can be rewritten into a more convenient form 



E=^J d 2 r J & 2 r'p(r)v(r, r')p(r') - -Nv(0) 



(7) 



where the integrations over r and r 1 are now taken over the whole space. 

The thermodynamic characteristics of the system are determined by the grand parti- 
tion function 5, 



oo oo N+ N- 

/6_l_ /J 



N + =ON-=0 N + ] - N - ] 



Q(N+,N. 



with 



r N 

Q(N + ,N_) = J IJdVxp [-pE({ qj , rj })} 

3=1 



3a) 



(8b) 



being the canonical partition function of N + positive and iV_ negative charges and iV = 
N + + iV_. To express S in terms of a 2D Euclidean sine-Gordon theory, we first recall 
that — A/(27r) is the inverse operator of v(r) [see eq. (0)]. The standard identity then 
follows 

/ £>0exp [/ d 2 r (f 0A0 + iy/W]3(f>p 



cxp 



d r / d r'p(r)v(r, r / )p(r J ) 



(9) 



/ V(f)exp ( J d 2 r| 

where <p(r) is a real scalar field and / Dip denotes the functional integration over this field. 
Inserting p from (|6]), the second term in the action of the field theory (P) takes a nonlocal 
form Qj[4>( x j, Uj) ~ 4 > (~ x jjVj)]- ft is therefore convenient to reformulate the field 

theory (Q) as a boundary problem using a procedure proposed in ref . [0 . One introduces 
two new fields 



V2 



(x,y)+(f)(-x,y)) 



{x,y) = -j= [(f>(x, y) - (f)(-x, y)] 



(10a) 
(10b) 



defined only in the positive x > half-space. Clearly, the even field has Neumann bound- 
ary conditions d(j) e (x , y) / dx\ x= o = and the odd field has Dirichlet boundary conditions 



b (x = 0,y) = 0. It is straightforward to show that 
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/ d 2 r^A0 = - / d 2 r {<j) e A<j) e + <p A<j) ) 
J 2 2 Jx>o 



11) 



2 r r 2 

The even field contributes to the action only by its free-field part e A0 e /2, which is 
"cancelled" with its counterpart in the denominator of (pj). Thus, when <p is renamed as 
</>, the rhs of (01) is expressible as a local field theory formulated in the half-space x > 0: 



/X?0exp{/ I>o d 2 r 




|(V0) 2 + V27r/3£ i <&0(r i )_ 


} 


JV(f)exp 


-4> o d 2 r|(V0) 2 





(12) 

with Dirichlet boundary conditions <p(x = 0, y) — 0. Now, defining by (• • •) the average 
over the field theory (|12"D , one proceeds along the standard line, i.e., express H as follows 

N+ / r . — \ jv_ \ 

: — — — ((/ ■■; 

JV + =OiV_=0 

(13) 



oo oo ^N + +N-J3ip(N + -N-) 
i—i v — > v — ^ ^ ^ 



where z = zexp[/3t>(0)/2] is the fugacity renormalized by the self-energy term, and after- 
wards sum over N, and N_. The final result reads 



JV(f)exp{f x>0 d 2 r 


-|(V0) 2 + 2^cos(v/2 


irf3(f) - i(3(p) 


} 


exp 


-L>od 2 ri(V0) 2 ; 





(14) 

with the fixed value of the field at the boundary, <p(x — 0, y) — 0. 

The field theory in (II) can be viewed as the ordinary 2D Euclidean sine-Gordon 



model in the half-space, defined by the action 

1 



A 



sG 



x>0 



(V0) 2 + 2^cos(/3 s ^ 



with 



(3 sG = ferti 

and Dirichlet boundary conditions <fi(x — 0, y) — (f>o, 

T 



(15a) 



(15b) 



2lX 



(15c) 



The boundary value 0o usually enters into the formalism in the combination 

e= ^-0o = -2i<p (16) 

PsG 

which will also be used in what follows. For if real (which is the case of interest), we 
assume that the solution of the present theory corresponds to an analytical continuation 
of the results of the theory with real Dirichlet boundary conditions. 

Without going into details we mention that one can proceed as above also in the case 
of the 2D TCP in contact with a dielectric of dielectric constant = 0. In that case, the 
images have the same charges as the source particles, the even field (|10a|) becomes the 



relevant one, and that gives a mapping onto an equivalent boundary sine-Gordon model 
with Neumann boundary conditions. 

3 Boundary sine-Gordon model 

The integrability property of the bulk (infinite in both x, y-directions) sine-Gordon model 
along the standard lines of the Bethe ansatz technique is well known from the seventies || . 
The discrete symmetry of the theory — > + 2im/ (3 s g {n integer) is spontaneously broken 
in the domain < 0^ G < 87r; one has to consider one of infinitely many ground states 
{|0 n )} characterized by (<f)) n = 27m/ (3 S G, say the one with n — 0. The spectrum of particles 
involves solitons S, antisolitons S and soliton-antisoliton bound states (breathers) {Bj,j = 
1,2,... < l/q} with masses 

raj = 2M sin (^if) (17) 



where M is the soliton mass. The parameter q, defined by 

8tt - [j sG 

ranges from to oo: the breathers exist for q G (0, 1) (or in the plasma stability range 
< (3 < 2) and there are no breathers for q larger than 1 (2 < (3 < 4). 



In the seminal work || it was argued that the restriction of the sine-Gordon model 
to the half space (x > 0, y) does not break its integrability if one adds a boundary action 
term 



A 



G , 



>B 



(19) 



(A s g)b = dy mcos 

J — oo 

where m and 0o are f ree parameters, and ^b{v) — <P{x, y) \ x =o- The underlying sine-Gordon 
theory (|T5|) corresponds to the m — > oo limit of (0) which sets 0b to 0o- 

The free energy of the theory ([15|) is related to the ground state energy of the boundary 
quantum (1 + l)-dimensional sine-Gordon model, defined by Lagrangian 

"I 



sG 



dx 



-(dcpy + 2zcos(#. 



sGS 



(20) 



L — ► oo, with boundary conditions 0(0) = 0o and = (fi' considered in terms of £ 
and respectively [see eq. (0)]. There exists a lattice regularization of the theory (|2H[), 



namely the XXZ model in boundary magnetic fields ]IU|, defined by Hamiltonian 



JV-l 



27T sin r 



[ a j a j+i + "i^j+i ~ ( cosr ) cr j cr j : +i] + hoi + h'o 



(21) 



with r G (0, tt/2). Here, e = 1 and e = —1 correspond to the ferromagnetic and antiferro- 
magnetic cases of the XXZ-chain, respectively. This model results as the hamiltonian limit 
of the inhomogeneous 6- vertex model on an open strip pT| , with an alternating imaginary 
part ±iA added to the spectral parameter on alternating vertices. The continuum scaling 
limit is given by taking A — >• oo, N —>■ oo, and the lattice spacing a — > 0, such that 
L = Na remains finite. In the bulk, the regularization fixes \JH oc (1/a) exp(— const x A) 
and 



sG 



sG 



8r (e = 1) 
8(tt - r) (e 



(22a) 
(22b) 



As r G (0, vr/2), the ferromagnetic regime corresponds to q G (0, 1) and the antiferromag- 
netic regime to q G (1, oo). As concerns the interrelation between the "surface" quantities 
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I2fl, one defines the function 



and 



T 

Then, introducing an auxiliary variable 



, /sinh i6-a/2\ . , 

f(a, b) = -i In — ^4 23a 

JK ' Vs nh 6 + a /2 J v ; 



H = -/(ir, -i ln(/i + cos r)) (23b) 



r 

one has 



t = - t>2 (24) 



e = \{t-l-H) (6=1) (25a) 



7T 



H 



« = 2 1 1 " J ( e = -D < 25b > 

The same formula hold for ^' in terms of if', resp. ft/. 

To study boundary effects in the XXZ model (0]), one looks for the solutions of the 



Bethe equations which correspond to a wave-function localized at j = or j = N and 
exponentially decreasing away from the boundary. These boundary bound states were 



identified with new boundary strings in the Bethe ansatz p4]| (for a review see 
In the thermodynamic limit iV — > oo, when the left and the right boundaries can be 
treated independently and the overlap of the corresponding wave-functions is negligibly 
small, the lattice ground state energy was found as a function of the parameter A of the 
related 6- vertex model. In the continuum limit A — > oo corresponding to the (1 + 1)- 
dimensional sine-Gordon model (|20|), the ground state energy of the last was obtained in 
the form 

Eground = -^feu/fc + Eyries + 0(1/ L) (26) 

where 



M 



2 



7T 



Ehuik = —L tan — (27a) 

4 \2(t-l)J v ; 



E> 



bdries 



M 

T 



sin (Hn/[2(t- 1)]) sin (H'n/[2(t - 1)]) 
sin(t7r/[2(t-l)]) + sin(*7r/[2(t- 1)]) 



cot 



for the ferromagnetic case e = 1 |T3| (note a different notation in this reference) and 



(27b) 



Ebulk 
Ebdries 



T M 2 Stir 
L—— cot — 
4 V 2 



M 

Y 



sin ((* - F)tt/2) sin((t- H')tt/2) ^ Stir 
sin (tn/2) + sin(tvr/2) C ° VT 



(28a) 
(28b) 



for the antiferromagnetic case e = — 1 [[H]]. In terms of the sine-Gordon parameters, 
"bulk" q (0) and "boundary" f,f both Eqs. (0) and Eqs. take the same 

form 



Ebulk 
Ebdries 



M\ fqn 

—L tan — 

4 V 2 

M [ 1 r f j-\ f -i i /?7T 

-y \c^(^) [cos(g0 + cos(ge } " 1 + tan It 



(29a) 
(29b) 



4 Surface tension of the plasma 

For a Coulomb gas of volume V bounded by a surface of area S, the grand potential 
Q = — /3 _1 lnS is the sum of a volume part and a surface part: 



n = -Vp(z,(3) + S 1 (z,f3) 



(30) 



where p is the pressure and 7 the surface tension. For a strip L x i?, R —>■ 00 and L large, 
the "specific" is given by 



lim " = -Lp(z,/3)+ 7 (^/5)U=o + 7(^/3)U=L + 0(l/L) 
The thermodynamics of (|i"4|) is mapped onto the ground state of (p0| ) according to 



(31) 



RE g rOUn d 



(32) 
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In the considered L — > oo limit, the boundary energy ( 29b ) is the sum of two clearly 
separated contributions comming from the boundaries at x = and x = L. To calculate 
the surface tension, one keeps only the contribution at x — 0, and identifies £ with the 
bulk potential cp via eq. fll6l). Thus, one gets 



Inserting ( [L5b ) into QTS[ ) and considering 
pressible in terms of the Coulomb ones as follows 



(33a) 
(33b) 

the sine-Gordon parameters are ex- 



Q 



Z = -2ip 



(34) 



4-/?' 

As concerns the link between the soliton mass M and the fugacity z, the formalism of 
section 2 showed us that z renormalizes multiplicatively. To give z a precise meaning, one 
has to fix the normalization of the field cos(/3 s g0)- The conformal normalization proposed 
in refs. [16[] and [|17] corresponds to the short-distance limit of the two-point correlation 
function 

(cos(& G 0)(z) cos(& G 0)(2/)> - \\x-y\-^a/^) (35) 

This normalization, equivalent to a well known leading short-distance behaviour of the 
positive-negative pair correlation in the Coulomb gas, fixes the z-M relationship as follows 



iWOz + i)) 



M 



V^T((?+ l)/2) 



1 2/(9+1) 



7rr(l/(g + l)) f" 2T(g/2) 
The total particle number density, generated via 

9^ 



(36) 



n = z- 



is related to M as follows 



n 



-M 2 {l + q) tan 



7TQ 



(37) 



(3? 
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The n-z relationship is given in paper I, Eqs. (49), (50). The singular behaviour of n as 
f3 — > 2 (z fixed) can be deduced from these formulae: 

4z 2 7T 



n ~ 



2-/3 



Finally, we have 



f3p = n I 1 



and 



n(4 - /?) 



2sin(7r/?/(4 - /?)) 



1/2 



+ sin 



2(4-/3), 



— cos 



Tf/3 



2(4-/3) 



/?7 has the small /3-expansion 



f3 1 = - (2vr/3n) 1/2 i 1 + 



7T 2<£ 2 



16 7T 



/3 + 



71" / 7T\ <£> 

64 V + 6 / + 2^ 



/3 2 + .. 



(39) 



(40) 



(41) 



(42) 



With regard to 
and 



for fixed, p and 7 exhibit the same type of collapse singularities as 



P 



Z 2 7T 



2-/3 



7 



2; COsh(2(y9) 

2-/3 



(43) 



(44) 



5 High-temperature expansion 



As a check of the exact expression fl41|) of the surface tension, the beginning of its high- 



temperature expansion [expansion in powers of (3, eq. fl42|)l will now be compared to a 
direct evaluation of the two first terms of this expansion derived from the renormalized 
Mayer expansion of the free energy of the Coulomb system. In the following, all functions 
and integrals are defined in the half-space x > 0. 
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The surface tension 7 can be defined as the boundary part per unit length of the 
grand potential Q. The total numbers of positive and negative particles, respectively, are 
N + = —[3z + dQ/dz + and N_ = — j3z_dVt/ dz_. Going to the variables z and ip defined by 
z± = ze ±/3lp gives for the total number of particles N = N + + iV_ = —(3zdQ/dz. The 
boundary part of this relation is 

— (3z— = dx \n(x) — n] (45) 
oz Jo 

Since, as recalled in I, z is proportional to n 1_/3 / 4 , ( [45]) can be rewritten as 



n ^7 

P n JT 
on 



1 ~a)Io° dX ~ ^ ^ 



This relation (^) will be used for computing the surface tension 7 from the density profile 
n(x) which will be determined as a function of the bulk density n, the electric potential 
(p in the bulk, and the inverse temperature f3. 

Our starting point is the relation obeyed by the density profiles: 

n ± (l) 



In 



z±(l) 



where A[n] is the negative of the excess free energy times /3, considered as a functional of 
the position-dependent densities, and z±(l) is the position-dependent fugacity at point 1. 
Since a particle at a distance x from the boundary has an interaction (l/2)ln(2x) with 
its own image, 

z± (x) = z exp [±/3<p - {(3/2) \n{2x)\ (48) 

The relation (^) is exact. Here, the renormalized Mayer expansion, described in I, is 
used for expanding A[n] up to order /3 2 , i.e., we keep only the contributions D (two field 
circles connected by a simple — f3v bond plus the sum of ring diagrams) and the graph 
D\. However the latter one will be shown to give no contribution to the density profiles at 
the desired order, and we shall be left with only D , which constitutes an approximation 

13 



of the Debye-Hiickel type. However, here, the interaction v(l,2) between particles 1 and 
2 includes the contributions from the images, i.e., v(l,2) = — lnri2 + lnr* 2 where 7*12 is 
the distance between points 1 and 2, r\ 2 is the distance between point 1 and the image 
of point 2. The contribution of Do to the functional derivative in fl4~7|) is 



=T(3 J d2 „(i, 2)[n + (2) - n_(2)] + ±[K(1, 1) + (3v(l, 1)] (49) 
where K(l, 2) is the renormalized bond defined by the integral equation 

K(l, 2) = -/3v(l, 2) + J d3[-/3v(l, 3)]n(3)if (3, 2) (50) 



This integral equation (|50| ) can be transformed into a partial differential equation by 
taking the Laplacian with respect to 1, Ai, of both sides of fl50|) , and using Aiu(l,2) = 
— 2tt8(1,2), which gives 

AxK(l, 2) = -2tt/M(1, 2) + 27r/3n(l)if(l, 2) (51) 

Since t>(l,2), the Coulomb interaction in presence of a conducting wall, vanishes when 1 
is on the wall, the same boundary condition holds for K(l, 2). 

The above equations can be solved for the density profile n(x) by iterations, starting 
with the lowest-order approximation of a constant n(x) = n in ([51]). Then, the solution 
of (|5T]), with its boundary condition, is 

K^(l,2) = -(3K (Kr 12 ) +/3K (Kr* 12 ) (52) 

where k 2 = 2ir/3n (k, is the inverse Debye length), and Kq is the modified Bessel function 
of second kind. Using this lowest-order K^> in (|4TS|), approximating A[n] by D [n] and 
using fl4lf ) in (|^]), gives 

n±(x) = ^exp (±/% - p(a;)] + 77 lim [-/^ (m~ 12 ) - /?lnr 12 ] + ^W (2kx)) (53) 
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where, with point 1 at a distance x from the boundary, 



d2v(l,2)[n + (2) -n_(2)] 



(54) 



(p(x) is the electric potential created at point 1 by the charge distribution (localized near 
the boundary) n + (2) — n_(2). ip — <p{x) is a finite quantity which goes to zero as x — > oo. 
Since, at the Debye-Hiickel order of approximation, the bulk fugacity and density are 
related by 

n = 2zexp{(l/2) lim \-(3K Q Ur 12 ) - p\nr 12 ]} (55) 
(j53|) can be rewritten as 



n 



n±(x) = -exp{±/% - <p(x)\ + ((3/2)K (2kx)} 



(56) 



(0) has a simple physical interpretation: each particle feels a mean one-body potential 
made of two parts: ip(x) is the electric potential created by the surface charge density 
in the plasma, —(1/2)Kq(2kx) is the screened interaction of the particle with its image. 



The linearized form of (|5q) gives for the density profile 



n(x) = n+(x) + n-(x) = n 
In the following, the integrals 



1 + -(3K (2kx) 



POO 

/ K (x)dx 
Jo 



71 



(57) 



(58a) 



and 



/•OO jr Z 

j o K 2 (x)dx = - (58b) 

will be needed. Using (|57| ) and (|58a|) in ( ^6|) gives the surface tension 7 at lowest order 
in (3: 



Pi 



: (2vr/3n 



,1/2 



(59) 



At this order in j3, an explicit form of <p(x) was not needed for computing the surface 
tension, neither does that surface tension depend on the parameter ip. However, (p(x) will 
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be needed in the following. It can be easily obtained by writing, from the linearized form 
of (|56|), n + (x) — = (3n[ip — <f(x)] and using the Poisson equation d 2 [<p — ip(x)]/dx 2 = 

2n[n + (x) — n,-(x)] with the boundary condition <p(0) = 0. One obtains 

ip — <f(x) = <£>exp(— kx) (60) 

The next iteration for n(x) is obtained by using the density (|56|) in the equation (|51|) 
for K and treating n(x) —n= (1 /2)nf3KQ(2nx) as a perturbation. Now K = + 
where is defined by (0). To first order in the density perturbation, ([51]) and (|52| ) 
give 

(Ai - K 2 )K m (l, 2) = 27T(35n(l)K^(l, 2) (61) 

where 5n(l) = (l/2)nf3K (2Kx) with x the distance of point 1 to the boundary. The 
solution of (|6l|), with the boundary condition that K^\l,2) vanishes when 1 is on the 
wall, is studied, by the method of Green functions, in the Appendix, where it is shown 
that K^Hl, 1) is a function of the coordinate x\ of 1 such that 

It is a priori necessary to keep also the contribution from the graph Di to A[n]. However, 
at the present order in /3, 5Di[n]/5n±(l) can be evaluated for constant densities and then 
it vanishes. Taking also K^ 1 ' into account, we now have, instead of ([56]), 

n ± (x) = | exp j±/% - <p{x)\ + ^K (2kx) + ^ (1) (^)]} (63) 

where K^'(l, 1), when 1 has the coordinate x, is renamed K^'(x). Expanding the expo- 
nential in (|63"D to order (3 2 gives for the total density n + + n_ 

n(x) -n = ^K (2kx) + \k^\x) + ^K^kx) + ^ - V (x)] 2 (64) 
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where it is sufficient to use for ip — <p(x) the lower-order expression (|60|). Thus, using (|58| 



and (|62T), one finds 



/ dx fn(x) — n] = — 
Jo 16 



1 /3 / 7T\ ^tt 2/3 2 
1 + -j I 1 — — ) + -3 — h — V? 
4 V 4/ 8 vr 



(65) 



Using (|65|) in (f46|) gives the final result 



/3 7 = _I( 2 7T/3n) 1/2 



(66) 



in agreement with (E2| 



6 Collapse singularity 

To check the singular behaviour of the surface tension 7 near the collapse point (3 = 2, 
eq. (|44T ), one starts with the exact x — > limits 

n+ ^- n+ ~j^m> n -( x) ~ n ~~j2^ (67) 

which can be derived directly by using the grand canonical or canonical formalisms, in 
analogy with the short- distance expansion of the positive-negative pair correlation in the 
bulk. At (5 = 2, the exact result for the density profile reads M 



OO 



n±(x) -n± = ^- I dl 

ATI JO 



m Ki exp(±/3<£>) + m 



exp(—2Kix) (68) 



Ki m cosh(/3<y?) + Ki 

where m = 2nz and K\ = (m 2 + Z 2 ) 1 / 2 (note that in the original work [|J there are some 
mistakes in the equation (3.25)). The x — > limit of eq. fl6"8"|) 



n±(x)-n±~^ (69) 
still is of the form (pTj), even when the n± densities diverge. Consequently, one can put 

nW-,= ^M/(2x) (70) 
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where / is a function regular in (3 around (3 = 2, with 



/(0) = 1 



(71a) 



The density n(x) is supposed to tend to its asymptotical x — > oo value n faster than any 
inverse power of x, so that 



lim f(2x) — > faster than any inverse power of 



(71b) 



According to (|45|), it holds 



c?7 
<9z 



zcosh(^) / m- p/2 f(t) 

Jo 



An integration per partes gives 



JKJ 1-/3/2 



du i-0/2d f(t) 
t=o J o dt 



(72) 



(73) 



For (3 < 2, t l PI 2 /(£) vanishes at t = as well as in the limit t —> oo due to the fast decay 
of f(t). Then, 

(74) 



^ _ cosh(/?y) [°° dttl ^ /2 df(n 



dz 1-/3/2 7o ^ 
When /3 — > 2~, one can perform a (2 — /3) expansion of the integral in (]7 



°° dt t 1 "^ 2 ^ = /(oo) - /(0) + 0(2 - 0) 



With regard to ( 71a|) and ( |71b|) , one arrives at the desired formula 



(75) 



7 Conclusion 



A two-dimensional model for the interface between an electrolyte and an electrode has 
been considered: The 2D TCP bounded by a rectilinear ideal conductor wall. Previously, 
the surface tension 7 in this model was known H only at the special inverse temperature 



(5 = 2 (in which case, for obtaining a finite result, a hard core repulsion between the 
particles and the wall had to be assumed). Now, the main result of the present paper, eq. 
(p|), provides the surface tension for point particles at any temperature, in the stability 
range of the model (3 < 2. 

The surface tension depends on the bulk density n as n 1 ^ 2 , as a priori expected for 
dimensional reasons. Its temperature dependence is more complicated than the one of the 
bulk pressure but simpler than the temperature dependence of some bulk thermodynamic 
quantities derived in I. 

Appendix 

In this Appendix, the correction to the renormalized bond is studied. In terms of 
the Green function which obeys 

(Ai - /t 2 )K (0) (l, 2) = 2*05(1, 2) (Al) 

with the boundary condition K^(l, 2) = when point 1 is at x = 0, the solution of (|6l|) 
with the same boundary condition is 

K^(l,2) = J d3ir (0) (l,3)K (0) (3,2)5n(3) (A2) 

where K^°> is given by fl52p and 5n(3) = (l/2)nl3K (2KX a ). Thus, the desired integral 

roo roo roo roo 3fl 

/ dx!^ (1) (l,l) = / dxt / dyj dx 3 [-/3K (Kr 13 )+(3K (Kr* 13 )} 2 ^-K (2Kx 3 ) 
Jo Jo J-oo Jo 2 

(A3) 

The integrals on x\ and y 3 will be performed first. Since the integrand is an even 
function of xi, the integral on x\ and y 3 is half that integral performed on the whole 



19 



plane. Furthermore, ys can be shifted into y\. Therefore, 

roc roo p 

I(x 3 ) = j o d Xl J ^dy 3 [-K (Kr 13 ) + K (Kr* 13 )} 2 = J d 2 n [K^r 13 ) - K (Kr 31 )K ( K r* 13 )\ 

(A4) 

One now has a convolution integral which can be performed by going to Fourier space: 

J(X3) = J dk w^W " 

= ^[\-2kx 3 K 1 {2kx 3 )} (A5) 

K 

Using ( |A5| ) in ( |A3| ) gives, using as a rescaled integration variable x = 2kx 3 , 

roo poo 

/ dx!A (1) (l,l) = / dx[l - xKAx)\Kq(x) (A6) 

JO oft JO 

After an integration by parts using K\{x) = —dK (x)/dx, and taking into account Q58b| ), 
one finds 
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dn 
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K- 
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